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We analyze the scaling properties of inclusive hadron production in proton-proton 
and in heavy ion collisions from fixed target to collider energies. At large trans- 
verse momentum pT, the invariant cross section exhibits a power- like behavior 
Ed^a/d^p oc p^"^ at fixed transverse x, xt = , and fixed center-of-mass 
scattering angle 9cm- Knowledge of the exponent n allows one to draw conclusions 
about the production mechanisms of hadrons, which are poorly known, even at 
high pt- We find that high-pr hadrons are produced by different mechanisms 
at fixed-target and collider energies. For pions, higher-twist subprocesses where 
the pion is produced directly dominate at fixed target energy, while leading-twist 
partonic scattering plus fragmentation is the most important mechanism at collider 
energies. High-pj' baryons on the other hand appear to be produced by higher-twist 
mechanisms at all available energies. The higher-twist mechanism of direct proton 
production can be verified experimentally by testing whether high pT protons are 
produced as single hadrons without accompanying secondaries. In addition, we find 
that medium-induced gluon radiation in heavy ion collisions can violate scaling. 
PACS: 13.85.Ni, 24.85.+p 



2 



I. INTRODUCTION 

A fundamental feature of quantum chromodynamics (QCD) and asymptotic freedom is 
the nearly scale-invariant behavior of quark and gluon two-particle hard-scattering processes. 
If these pointlike hard-scattering subprocesses are convoluted with the structure functions of 
the incident hadrons and the fragmentation functions which produce final state interactions, 
the resulting inclusive cross section E "^ (T{hahb-*hx) g^g^jgg nominally as at fixed xr = 
^^^f^ and 9cm (We denote the center-of-mass (cm.) scattering angle by 9cm and the three 
momentum of the produced hadron in the cm. frame by Pcm- The cm. energy is ^/i)- order 
to validate this scaling at 9cm = 90°, one needs measurements at various cm. energies and 
corresponding values of Pcm- The conformal scaling prediction is modified by the logarithmic 
running of the QCD coupling and the logarithmic corrections to the scale-free parton model 
arising from the evolution of the structure functions and fragmentation functions. These 
results are rigorous consequence of perturbative QCD factorization theorems for inclusive 
hadron reactions at large p^. 

In the same way as Bjorken scaling provides evidence for quarks in deep inelastic lepton 
scattering, the scaling behavior of the cross sections for the production of high particles 
in hadronic collisions can be used to test the scaling of the underlying QCD subprocesses as 
well as the existence of point-like constituents each carrying a finite fraction of the hadron's 
momentum . However, in some cases such as hard exclusive reactions, the underlying 
hard subprocess cannot be the simple 2 — 2 reactions; since all of the valence quarks of the 
interacting hadrons are involved in the scattering process. Even in inclusive reactions, higher 
parton number processes can contribute, particularly at high xr where there is a trade-off 
between the fall-off at high px which favors minimal number of scattering partons and the 
fall-off at (1 — xr) which favors more valence partons entering the hard subprocess. 

It is clearly important to carefully analyze the scaling behavior of inclusive reactions 
in order to confirm the validity and applicable kinematic range of the leading-twist 2^2 
subprocesses. Our aim is to determine the mass dimension of the microscopic matrix element 
from dimensional and spectator counting rules , which have recentl 

nonperturbatively for nearly-conformal theories using AdS / CFT duality 
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been derived 
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Since the dimension of this matrix element is related to the number of participating ele- 
mentary fields, such an analysis provides detailed information about the specific microscopic 
hard process underlying high pt hadron production. This is particularly important for the 
interpretation of hard processes in nuclear collisions. 



II. DIMENSIONAL COUNTING RULES 



active 



The partonic 5* matrix element is related to the partonic amplitude by Sfi 
Sfi + «(27r)^5*^^)(^pj„ — J2Pout)Afi- With single-particle states normalized to {p\p') 
2Ep(27r)^5^^\p — p'), the amplitude Afi on the microscopic level has dimension mass^ 
Here UacUve is the number of active fields, i.e. the number of elementary fields entering Afi. 
The simplest example is 2 ^ 2 scattering between quarks and gluons. In this case, UacUve = 4 
and the partonic matrix element is dimensionless, as is natural for a scale invariant theory. 
However, because of color confinement the partonic subprocess cannot be observed directly, 
and one needs to find a way to connect to the hadronic reaction. 

The kinematics of an inclusive reaction hahb hX is described by 3 Lorentz invariants 
These are e.g. the center-of-mass energy squared s = {Pa + Ph)'^, the transverse momentum 
transfer (squared) t = {Ph — PaY the missing mass Mx. It is common to introduce the 
dimensionless variables {u = M]^ — s — t), 

u t 

Xi = , X2 = — . (1) 

s s 

These variables are related to the rapidity y and radial xr of the observed hadron by 

y = hog(-), (2) 



2 °\X2 

XR = %^ = l-Ml/.^xi + X2, (3) 

Since most of existing data are at y = where xr = xt = 2pt/\^, one often refers to the 
scaling of the invariant cross section as "xy scaling". For y ^ 0, we find the variable xr 
more useful than xt, since xr allows a smooth matching of inclusive and exclusive reactions 
in the limit xr ^ 1. 
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We shall assume that at high px-, the inclusive cross section takes a factorized form, even 
if the microscopic mechanism is higher twist, 

da{hahb hX) = ^Ga/ha{xa)Gb/ht{xb)dXadxb^ \Afif dXfDh/c{zc)dzc. (4) 

abc 

The dimensionless functions Ga/hai^a) describe the momentum distributions of partons of 
type a in hadron ha, where a may stand for quarks and gluons as well as for composite degrees 
of freedom, such as diquarks and intrinsic hadrons. These functions cannot be calculated 
perturbatively, except in the limits Xa,b 1- For quarks and gluons, the scale dependence of 
the distribution functions is described by the DGLAP evolution equations, but the evolution 
of color-neutral degrees of freedom is suppressed by at least one power of the hard scale, 
since gluon radiation off color neutral objects is suppressed. Similar observations can be 
made for the fragmentation function Dh/c{zc), which accounts for the transition of a parton 
c into a hadron h with momentum fraction Zc = Xi/xa + X2/xb- The amplitude of the hard 
subprocess Afi is assumed to be calculable in perturbative QCD. Integration and summation 
over all unobserved variables, such as the phase space dXf of the final state, is understood. 

By keeping all ratios of Mandelstam variables fixed, the x dependence of the distribution 
functions does not affect the scaling behavior of the hadronic cross section. The factorization 
hypothesis Eq. then yields the power law 



JMhghb ^ hX) f{t/s,u/s) 
d^p 



^ yU. = n -2 ' (5) 

w, »,v^ e ' ''active ^ 



which reflects the mass dimension of the microscopic amplitude. Hence, the inclusive cross 
section multiplied by with 

n = 2nactive - 4, (6) 
is a function of the dimensionless variables y and xpt only, 

^ d^ajhah hX) ^ F{y,XR) 
d^p p^^^'^-^^ ' 

This is the desired relation: the pr dependence of the inclusive cross section is directly related 
to the number of participants UacUve in the microscopic matrix element. In higher twist 
processes, the function F{y, xr) also depends on the hadron distribution amplitudes, which 
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have dimension mass for mesons and dimension mass squared for baryons. For example, 
the pion distribution amphtude is normahzed to 7,^+ = 130 MeV. The large-p^ behavior of 
Eq. ((7j) is determined by the minimum number of active partons associated with the reaction 
ha + hb h + X. Of course, n will depend on y and xr- Let us illustrate the scaling behavior 
by the following two examples: 

• For pp — i> pX, one can have the subprocess uu —>■ uudd, which has riacuve = 6 rather 
than 4 active elementary fields. Because of proton compositeness, the inclusive cross 
section E^-^^^^^^^ scales nominally as ^ at fixed xr and y where the dimensional 
factor /tv reflects the physics of the proton distribution amplitude. In this process, 
the proton is made directly in the short distance reaction rather than from quark 
fragmentation or resonance decay. Because no energy is wasted in the fragmentation 
process, the cross section falls off relatively slowly ~ (1 — xj^y at large xr. In general, 
one finds the power law behavior (1 — x/j)^'"^"^ of the inclusive cross section at large 
Xr, where is the total minimum number of spectators in ha, and h, which do not 
participate in the hard scattering reaction. In this example = 2 + 2 = 4. [We ignore 
here complications from parton spin.] 

• In the limit xr — > x^'^^ ~ 1, the missing mass Mx reaches its minimum allowed value 
and the reaction becomes exclusive. In this case, there are no spectator fields, so 
that the number of active participants UacUve attains its maximum value, e.g. n = 
'inactive — 4 = 20 for pp pp. Thus, the leading-twist pQCD approach to high pt 
hadron production must fail as xr increases toward unity and the process becomes 
exclusive. From the definition of xr in Eq. Q it becomes clear that this is the case at 
large Xi and small X2 (or vice versa), i.e. for high pj- hadron production at large \y\. 

Finally, we remark that Xi^2 are not momentum fractions. The latter are denoted by Xa,b- 
The momentum fractions are defined only within a given parton model of hadron production, 
and cannot be related to hadronic invariants. Instead, the factorization ansatz for the cross 
section Eq. leads to integrals over Xa,b- In the case of 2 — 2 hard scattering, the lower 
integration limits are 
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1 — Xi 
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In the exclusive limit Xi^2 (1 i Xp)/2 [xp is Feynman x), both momentum fraction ap- 
proach unity, so that only valence partons are important. Hence, inclusive hadron production 
at very large rapidity is unaffected by gluon saturation. Such coherence effects disappear at 
the largest xp. (The authors of Ref. j8[ come to a similar conclusion from a different view- 
point.) This shows that the high energy limit of QCD cannot be completely described by the 
color glass condensate j^. However, (nearly) exclusive reactions at Aqcd Pt ^ still 
allow one to study perturbative QCD processes in a kinematic regime where Regge theory 
applies. 



III. PHENOMENOLOGICAL APPLICATIONS 



In real QCD the nominal power laws discussed in the previous section receive corrections 
from the breaking of scale invariance in QCD, i.e. from the running coupling and the scale 
breaking of structure functions and fragmentation functions. These corrections have been 
discussed a long time ago in Ref. [l^ but have not yet been studied quantitatively. 

Including scaling violations, the inclusive cross section of Eq. (0) changes to 

, 2n» / T 2 



d?a{hahb hX) 
hi — 



'^'^^ ^ ""^j^, c^T^ikijm- (9) 



Pt 



The threshold behavior of the cross section follows from spectator counting rules [10|]. We 
ignore here an extra contribution to this power which arises from helicity mismatch in the 
fragmentation process. The strong coupling constant a^*^" (^x^) (^s ^^e number of spectator 
fields) arises at large momentum fraction, since all spectators must combine their momentum 
to produce one high-x quark. This quark is far off-shell with virtuality = — so 
that the high-x tail of the structure function is calculable in perturbative QCD. Here, krjs 
the transverse momentum of the quark and fhg is related to the quark mass, see Ref. [l(j| 



for details. 

Eq. matches smoothly onto the exclusive limit Xr 1. This is still true in the presence 
of scaling violations: the correction to the simple power 2ns ~ 1 due to gluon radiation is 
contained in the function 

aPT) = — / ^c^sikl) = ^ In , ^J'^,y (10) 



7 



Heff 

6f 



5 . 75 
5 . 5 
5.25 

5 

4 .75 
4 . 5 
4.25 



10 



15 



20 



Pt (GeV) 



FIG. 1: The effective power n-eff according to Eq. Ml]) . The lower curve assumes 4 active fields 
and asymptotically approaches neff{nactive = ^^PT oo) = 4. Calculations were performed at 
XT = 0.03 and y = 0, which are typical values for RHIC. 

Here, (3q = 11 — 2Nf/3 is the QCD /5-function, Cr = Cp = 4/3 for quarks and Cr = 
Ca = 3 for gluons. Note the lower integration hmit k1^-. at large xr, the phase space for 
gluon radiation vanishes and QCD scaling violations disappear. Hence, the simple spectator 
counting rules become exact at the exclusive boundary. 

We shall now investigate, how QCD scaling violations affect xr scaling. For that purpose, 
we define an effective power Uefflpr) by taking the logarithmic derivative 

-"'^^^ 

of the cross section. 

We first concentrate on RHIC kinematics at ?/ = 0, where rather low values of xt ~ 0.03 
can be reached. Therefore, we drop all factors describing the large xr behavior of the cross 
section and determine rie// from the running coupling only. Different choices of the hard 
scale change numerical results by only few percent. We also checked that the (1 — xr) term 
is numerically irrelevant. Results are shown in Fig. ^ For the lowest order process 2 — 2 
process we find that the effective power rig// approximately increases by unity. This is close 
to what is seen in direct photon production at RHIC {ueff ~ 5) 

Following the suggestion of one of us (SJB), the PHENIX collaboration has analyzed 



the scaling properties of data 



12j |. For neutral pions, a value of n^ff = 6.33 ± 0.54 has 
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FIG. 2: Invariant cross sections for pp {tt^ + Tr~)/2 + X at three different energies (y/s = 19.4 
GeV, 23.8 GeV and 27.4 GeV) multiplied by pj,. The power nf.ff = 8 indicates a higher-twist 
mechanism. The curve shows the (1 — x/j)^ threshold behavior. 

been determined, which is somewhat larger than the leading-twist values shown in Fig. CI 



Strictly speaking, the analysis of Ref. 



12l | was done for the invariant hadron yields rather 



than the cross sections. Taking into account the variation of the inelastic cross section 
between -^i = 130 GeV and 200 GeV, the value of rie// turns out to be about 6.16, which 
is still within error bars. However, since next-to-leading (NLO) order perturbative QCD is 



able to reproduce RH 



C data on pion production, using fragmentation functions from e^e 



3 



we conclude that pion production is dominated by leading-twist 



annihilation as input 
processes. 

On the other hand, the Chicago-Princeton data exhibit a strikingly different power law 



1J|. An analysis of early 



and could never be described in the conventional parton model 
data on inclusive tt"*^ production yields n = 8.2 ± 0.5 for xr = xt > 0.35, i.e. pt ^ 3.5 



GeV at -^i ^ 20 GeV Similar results are obtained for vr , see Fig. |21 The power 
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FIG. 3: Protons produced in AuAu collisions at RHIC do not exhibit clear scaling properties in the 
available px range. Shown are data for central (t) — 5% ) and for peripheral (QQ — 90% ) collisions. 

law Ed^a/d^p{pp — » tt+X) oc p^^'^ giving Uactive = 6 may indicate a quark-quark scattering 
process which produces in addition to the incoming quarks a qq pair, which becomes the 
observed pion with high transverse momentum. This process has been analyzed within the 
Constituent Interchange Model (CIM) Q], where an incoming qq pair collides with a quark 
by interchanging a quark and antiquark. The CIM is motivated by the inclusive to exclusive 
iransition mentioned above and is in good agreement with the Chicago-Princeton (CP) data 
isl ] . The model even can reproduce the absolute normalization of the inclusive cross section. 

Obviously, the production mechanism for high px hadrons changes from -^i = 20 GeV to 
a/s = 200 GeV. For constituent interchange longitudinal momenta of 0(1 GeV) can still be 
accommodated in the wave function of the proton. When the relevant longitudinal momenta 
are about 0(10 GeV) at higher energies, interchange is no longer possible which the different 
reaction mechanisms with increasing energy. 

Moreover, for proton production the pt dependence at Chicago-Princeton energies is 
also explained by CIM. A value of n = 12 is a strong indication that higher twists from 
wave function effects dominate high pt hadron production around ^/s = 20 GeV. Here the 
produced proton is the result of proton scattering on a quark. If protons and pions were 
both produced by fragmentation as in the Feynman-Field-Fox parton model, it is hard to 
understand how a dimensionless fragmentation function could change n from 8 for pions to 
12 for protons. 
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Since high-pT- protons are produced by higher-twist mechanisms at fixed target energies, 
we also investigate the scahng properties of proton production at RHIC. The points in 



m and the 200 GeV data of Ref. 



Fig. 121 were obtained from the 130 GeV data of Ref. 
Unfortunately, the data do not extend out to large enough px and error bars become too large 
at high pt to establish xt scaling. It is important to measure inclusive proton production 
out to larger pt for at least two values of ^/s. From these data one could find out whether 
proton production is leading or higher twist. If protons are produced in nuclear collisions 
by parton recombination (see e.g. Q]), the cross section should fall off exponentially, i.e. 
there would be no xt scaling. 

IV. NUCLEAR EFFECTS 

It is interesting to investigate nuclear effects on the observed scaling laws, i. e to compare 
the scaling properties of Ed^a / dP'p{pp — > hX) to Ed^a / d^p{AA hX). Since pions appear 
to be produced by a leading-twist mechanism, the quenching of pion spectra may be due to 
medium induced gluon radiation. In the following, we shall adopt the formalism of Baier et 
al. (BDMPS-Z formalism, see Ref. for a review). 

The presence of a new dimensionful scale in nuclear collisions, namely the BDMPS trans- 
port coefficient q, gives rise to the possibility that xt scaling is modified or violated. With 
the quenching weight Q{pt) defined according to 

1 d^a{AB hX) d^a{pp -> hX) 



Ncoii dydp^ dydp% 
medium effects may modify rieff as 



-Q(Pt), (12) 



n-' = - ^^^^ (13) 
dmpT 



The calculation of Q{pt) has been performed within the BDMPS-Z formalism in Ref. 
These are the basic steps: Let P{AE) be the probability that a fast parton loses energy 
AE due to gluon radiation. The medium modified cross section can then be written as the 
convolution 

-^ = y4A^)P(AE)-^(p, + AE). (14) 
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Since the partonic cross section is a steeply falling function of px^ a small value of AE' 
produces a large suppression. For AE <^ pj-, one can expand the logarithm of da'""-'^ / dp\{j)T+ 
AE) in Eq. (HD) and obtains 3 

Q{Pt) ^ J d{AE)P{AE) exp (--—n^^'j = exp J duj{l - e"^)^ 1 • (15) 

The advantage of this approximation is that one obtains a particularly simple relation 
between Q{pt) and the gluon multiplicity dl/du, provided one also assumes a Poisso- 
nian probability distribution P{AE). The upper integration limit in Eq. (fT3j) is given by 
^^max = niin(a'LPM, -E'), where E = pt is the energy of the fast parton. Omitting overall 
numerical constants, the Landau-Pomeranchuk effect in QCD has the gluon spectrum 

— oca fl6) 

where ojlpm = {L is the length of the traversed medium) will be estimated below. 

Parametrically, we obtain the following result: depending on whether the energy p^ of 
the fast parton exceeds the critical value ujlpm, one distinguishes the two regimes, 



^ dlnQjpr) ^ ^su"-^^ PT « U.LPM ^^^^ 



dlnpj 



ttsn™'-™ Pt > UJLPM. 



Hence, xt scaling should be strongly violated in nuclear collisions, in contradiction to what 
is seen in experiment. Data are available from the PHENIX collaboration at RHIC, see 
Fig. 18 of Ref. Remarkably, despite the strong nuclear suppression of pion spectra, Ue/f 
for neutral pions has almost no centrality dependence, i.e. rieff = 6.33 ± 0.54 for peripheral 
and Ueff = 6.41 ± 0.55 for central collisions. We conclude that radiation of a large number 
of soft gluons is not the dominant mechanism behind jet quenching. 

Once again, we stress the importance of studying the cross section at fixed xt and ra- 
pidity, rather than at fixed ^/s. In the latter case, the pt dependence of the inclusive cross 
section is affected by the x-dependence of the structure functions, which can result in a 



Pr-independent quenching ratio 



21| . An analysis at fixed xt and y eliminates the sensitivity 



to parameterizations of structure and fragmentation functions. 

However, including all charged hadrons, the power n increases with centrality from n = 
6.12 ± 0.49 in pp -> + h^)X to n = 7.53 ± 0.44 in AuAu {h+ + h")X. The reasons 
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for this difference may be in baryon production, since the inclusive baryon cross section has 
a steeper pj- dependence. In heavy ion colhsions, pions are strongly suppressed while proton 
production at not too large pt has almost no centrality dependence. We therefore argue 
that n^ff = 7.5 reflects the scaling behavior of baryon production at RHIC. The underlying 
mechanism could be the process uu pd with Ueff = 8 as explained above. 

Finally, in order to estimate ojlpm (or equivalently g), we shall rely only on data that 
are not related to jet quenching. Using Bjorken's estimate of the initial energy density, one 
obtains 

eB, = ^ i'^-f] ^ 10 GeV/ fm^ ^ 606,,,, (18) 
at initial time tq = 0.5 fm. We account for the longitudinal expansion of the medium by 



employing the dynamical scaling law of Salgado and Wiedemann 
expanding medium to an equivalent static scenario, 

.hot 2g(ro) r+^^^,^ ^,,,eo^. ,GeV 



which relates the 



= -jf- J dr{r - To)f ^ lOr'" ^ (19) 

For L = 5 fm we obtain ujlpm = 25 GeV. The value of g™'"^ has been estimated in j^]. 
Even though our estimate of wlpm may be uncertain by a factor 3 in each direction, we are 
quite sure that RHIC high px data lie in the regime where the BDMPS-Z medium induced 
gluon radiation would yield an energy loss dE/dz oc —agy/qE, similar to the QED Landau- 
Pomeranchuk effect. A mean energy loss proportional to the projectile energy may be more 



consistent with the data 



24|. 



V. SUMMARY 



We have reviewed dimensional counting rules at xt scaling laws for high pr inclusive 
hadron production At leading twist, the inclusive cross section Ed^a/d^p scales 

nominally as p^" with = 4 at fixed ratios of invariants. Scaling violations in QCD, in 
particular the running coupling constant slightly increase the value of n. The experimental 
value rieff = 6.33 ± 0.54 for neutral pion production in peripheral heavy ion collisions at 
RHIC is somewhat larger than the expectation from leading twist. Nevertheless, we think 
that leading- twist partonic scattering is dominant, since NLO perturbation theory describes 
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the data reasonably well At fixed target energies however {^/s ^ 20 GeV), an effective 
power of rieff ~ 8 for pions and rieff ~ 12 for protons is a strong indication for higher- 
twist mechanisms in the SPS energy range. The fixed target data can be reproduced in the 
Constituent Interchange Model of Ref. Jj. 

Measurements of single particle hadron and photon production at the new hadron facilities 
at GSI and J-PARC will be very sensitive to higher twist effects. We have shown how 
one bridges the high and low energy domains. In addition, our analysis can be used to 
properly obtain the exclusive-inclusive connection. The conventional approach overestimates 
evolution at the exclusive limit. 

We have also investigated how medium induced gluon radiation changes the scaling prop- 
erties of high-pT- hadron production within the BDMPS-Z formalism. We find that radiation 
of an infinite number of soft gluons would violate Xt scaling. This is however not what is seen 
in experiment: the scaling law of neutral pions is unaffected by the nuclear medium (within 
error bars) [l^]- We conclude that medium-induced gluon radiation is not the mechanism 
responsible for pion quenching at RHIC. This is further supporte d by the fact that charm 



1 by tJ 
23, ^. 



e. t 



lere is no dead 
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production appears to be strongly suppressed in nuclear collisions 
cone effect as one would expect if quenching were due to bremsstrahlung 

Including all charged hadrons, the effective power law changes from rieff = 6.12 ± 0.49 in 
peripheral colhsion (60-80%)to Uejf = 7.53 ±0.44 in central collisions (0-10%) at RHIC. We 
argue that the larger value of Ueff in central collisions is due to a large baryon contribution 
to the charged hadron yield. If protons are produced by a higher-twist mechanism such as 
uu pd at not too high px, then there would be an intermediate px range in which protons 
scale with Upff ~ 8. This higher-twist mechanism is different from parton recombination 

n 

models [Iq], which lead to an exponential fall off. Unfortunately, existing data on inclusive 
proton production at RHIC do not allow one to draw definite conclusions at this time. Data 
at different energies will be required to separate the contributing mechanisms. In addition, 
the higher-twist mechanism of direct proton production which we propose can be verified 
experimentally by testing whether high pj- protons are produced as single hadrons without 
accompanying secondaries. It is clearly essential for the correct interpretation of the heavy- 
ion collision data that the role of the higher-twist processes be definitively determined. 
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